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We consider the entanglement entropy for a line segment in the system of noninteracting one- 
dimensional fermions at zero temperature. In the limit of a large segment length L, the leading 
asymptotic behavior of this entropy is known to be logarithmic in L. We study finite-size corrections 
to this asymptotic behavior. Based on an earlier conjecture of the asymptotic expansion for full 
counting statistics in the same system, we derive a full asymptotic expansion for the von Neumann 
entropy and obtain first several corrections for the Renyi entropies. Our corrections for the Renyi 
entropies reproduce earlier results. We also discuss the entanglement spectrum in this problem in 
terms of single-particle occupation numbers. 



I. INTRODUCTION 

Entanglement is one of the central concepts of modern 
quantum mechanics and quantum information theory. It 
characterizes the amount of correlations between parts 
of a quantum system. In recent years, a progress has 
been achieved in studying entanglement for a variety of 
models, with the most detailed results available for one- 
dimensional systems, see e.g. the review^. 

Entanglement can be introduced in a particularly sim- 
ple way in the case of a many-body system in a pure 
state, e.g., in the zero-temperature ground state, which 
will always be assumed in this paper. Let such a system 
be divided into two subsystems A and B. Then the en- 
tanglement may be characterized by the properties of the 
reduced density matrix pa of the subsystem A, which is 
obtained by tracing out the remaining degrees of freedom 



p A = tr B P 



(1) 



(here p denotes the density matrix of the pure state of 
the total system). The (von Neumann) entanglement 
entropy is then defined as the von Neumann entropy of 
Pa-, 



S {A) = -txp A hxp A 



(2) 



Though characterizing entanglement by a single number 
is appealing, it falls short in representing its full complex- 
ity. A more complete description of entanglement may 
be given by the set of Renyi entropies 



S (A *> = 



1 



1 - a 



In tr p^ , a > , a^l 



(3) 



(the von Neumann entropy can then be expressed as the 
limit S (A) = lim a ^iSi A) ). 

Since the total system is assumed to be in a pure state, 
these definitions can be shown to be symmetric with re- 
spect to the interchange of the subsystems A and B: 

S {A) = S (B) for both 

von Neumann and Renyi entropies^, 
so we shall drop the superscript (A) or (B) in our nota- 
tion below. 



Equivalently, entanglement may be characterized by 
the spectrum of the reduced density matrix pa (which 
coincides with the spectrum of ps for a pure statejP^. 
Like the full knowledge of the Renyi entropies, the en- 
tanglement spectrum allows to determine the state of the 
system up to unitary transformations in the subsystems 
A and B. In this sense, the Renyi entropies and the en- 
tanglement spectrum contain the full information about 
entanglement. 

The problem of calculating the entropies or the entan- 
glement spectrum simplifies in the case of noninteracting 
particles (bosons or fermions). In this case, the reduced 
density matrix {pa or ps) can be factorized into density 
matrices of individual single-particle levels^, and both 
the entanglement spectrum and the entropies may be ex- 
pressed in terms of single-particle occupation numbers Pi. 
In the case of noninteracting fermions, the entropies are 
given by 



S = - ^ \pj In ft + (1 - Pi) ln(l - pi 

i 

for the von Neumann entropy and 

S a = j^- ^ln[p? + (l- Pi r] 
1 — a — ' 



(4) 



(5) 



for the Renyi entropies. The sums over i can be converted 
into integrals over p, t [Eqs. (17) and (18) below] by intro- 
ducing the spectral density of the occupation number 



(6) 



This spectral density, together with the entanglement 
entropies Q and ([5]), in the model of noninteracting 
one-dimensional fermions, will be the main object of our 
study. 

Note that, in the case of noninteracting particles, the 
same spectrum of occupation numbers pi defines the full 
counting statistics (FCS) of the number of particles in 
each of the two subsystems. This observation was used 
in Refs. [7HH] to establish an exact relation between the 



2 



FCS and the entanglement spectrum. In the case of non- 
interacting fermions, both the FCS and the entanglement 
spectrum can be expressed in terms of the spectrum of a 
single-particle correlation matrix (in the context of FCS, 
such a decomposition was done in Ref. 1101 on the basis of 
the Levitov-Lesovik determinant formula^). 

Furthermore, for nonintcracting fermionic systems 
with translational invariance, the corresponding spec- 
tral problem involves matrices of Toeplitz type. There- 
fore, the asymptotic behavior of FCS and entanglement 
spectrum in the limit of a large subsystem size may be 
obtained with the help of the theory of Toeplitz de- 
terminants. A prominent example is the spin- 1/2 XX 
chairP^EOS which can be mapped to a system of noninter- 
acting fermions via a Jordan- Wigner transformation. In 
many interesting one-dimensional situations (including 
free fermions), the relevant Toeplitz matrix has Fisher- 
Hartwig singularities, and the asymptotic behavior of 
its deter minant can be found using the Fisher-Hartwig 
conjectur e! 14 * 15 *. While the leading asymptotic behavior 
of entanglement and FCS can be obtained by choosing 
the main Fisher-Hartwig branch, finding subleading con- 
tributions requires more work. Recently, corrections to 
the entanglement entropies accurate to order 0(L~ 3 ) (for 
a block of size L) have been computed for the spin-1/2 
XX chair{23 and in the continuous limilP^ 

Furthermore, in the continuous limit, a full asymp- 
totic expansion of the corresponding Toeplitz determi- 
nant was conjectured in Ref. [T5] in the context of FCS. 
Based on the matrix Riemann-Hilbert problem^ and, 
independe ntly, o n the Painleve V equation in the Jimbo- 
Miwa for m 1 20 ! 21 -*, an expansion was constructed for the 
FCS generating function of the particle number on a line 
interval for one-dimensional free fermions in the zero- 
temperature ground state. Using the periodicity conjec- 
ture for the expansion (not proven, but verified up to 
high orders in l/L), the asymptotic expansion was writ- 
ten in an explicitly periodic Fisher-Hartwig forrrP^l In- 
stead of selecting the leading Fisher-Hartwig branch, all 
the branches were combined to obtain a full expansion 
to all orders in l/L, taking into account the switching of 
branches intrinsically. 

We use the relation between FCS and entanglement 
entropies to carry over the full expansion conjectured in 
Ref. HU of the FCS generating function to the problem 
of finding the entanglement entropies and the entangle- 
ment spectrum for free fermions on a line. In partic- 
ular, we find the power-law asymptotic expansion for 
the von Neumann entropy S, compute first several co- 
efficients, and present an algorithm for calculating the 
coefficients to an arbitrary order. A similar approach to 
the Renyi entropies S a produces an expansion with oscil- 
lating terms. For the Renyi entropies, we only compute 
the lowest-orde r term s, which agree with the previously 
available results^EH. We also find finite-size corrections 
to the spectral density of single-particle occupation num- 
bers n{p). 

The remaining parts of the paper are structured as 
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FIG. 1. The function r(p) in Eq. pi. At the end points of 
the interval [0, 1], r(p) tends to infinity (logarithmically). 



follows. Section |TT] embodies our main results. In sec- 
tion |III[ we review the asymptotic expansion of the FCS 
for one-dimensional free fermions. The calculations of 
the spectral density fj,(p) and of the von Neumann and 
Renyi entropies can be found in section |rv] Our results 
are further discussed in Section[V] The appendix includes 
details of the analysis of oscillating terms in the asymp- 
totic expansions of the entanglement entropies. 



II. RESULTS 

Based on the conjecture for the FCS in Ref. [T5], we 
derive the asymptotic power series for the entanglement 
entropy of free fermions on a line in the ground state: 



S(x) = - ln(2x) + T + *2»s 



(7) 



71=1 



Here x = kpL (L is the length of the line segment for 
which the entanglement is computed and Uf is the Fermi 
wavevector) and the constant T is given by Eq. (31 1. 



Note that this series contains only even powers of TJx. 
All the coefficients S2n are rational numbers which can 
be computed to any given order in n using the methods 
of Ref. 18. The first several coefficients are: 



S2 



1 

12 



S4 



31 

96 



S6 



7057 
1440 



(8) 



The leading term (1/3) ln(2x), the constant t erm T, a nd 
the coefficient S2 are known from earlier works^EHSl. In 
contrast to the Renyi entropies, there are no oscillating 
contributions to the von Neumann entropy at any order 
in 1/x. 

The calculation involves an expansion for the spectral 
density fi(p) based on the conjecture in Ref. [TS1 Away 
from the end points p = and p = 1 (a precise condition 
is formulated in Section IV A), the spectral density has 



a quasiclassical structure with locally nearly equidistant 
levels. The smooth (nonoscillating) part of the spectral 
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density is given by 



fi{p, x) 
where 



1 



TV 2 p(l — p) 

t(p) = <p' 



[ln(2x) - r{p)] + 0{x- x ) , (9) 



2tt V p 



<P(0 = ar § 



(10) 



(11) 



and prime denotes the derivative of <^?(£) with respect to 
its argument. The function r(p) is plotted in Fig.[l] 



III. FCS OF FREE ONE-DIMENSIONAL 
FERMIONS 

We consider free spinless fcrmions on a continuous line. 
The temperature is assumed to be zero, i.e. the system is 
in the ground state characterized by the Fermi wavevec- 
tor kp. We will study the entanglement between two 
subsystems: an interval of length L and the remainder of 
the line. Both FCS and the entanglement in this setup 
depend only on the dimensionless parameter x — kpL. 
For example, the average number of particles on the line 
segment is given by (N) — x/ir. The FCS generating 
function for the probability distribution of the particle 
number N, 



J (27T«) N 



was conjectured in Ref. 1181 to be given by 

oo 

x(k,x) = ^2 X*{k + j,x) , 



(12) 



(13) 



X*(k,x) = exp 



2inx — 2k \iix + C(k) 



C{k) =ln [G(l + k) 2 G(1- k) 2 



2k 2 In 2 



(14) 
(15) 



where G(z) denotes the Barnes G- function and f n (n) 
are polynomials in k, computable order by order. For 
our purpose, we will use the logarithm of this expansion, 
which for— 1/2<k<1/2 takes the form: 



In x( K > x ) — 2?kx — 2k 2 In x + C(k) 



(16) 



OO |_ 2 J 

E E 

n=l m= _L ? j 



where |_ - J denotes the integer part of the argument. 

The coefficients C„ jm (/t) can be expressed in terms of 
the polynomials f n ( K ) order by order. They are also 
linearly related to the coefficients R n ,m(n) used in Ref.H8l 
for the expansion of the derivative (in x) of Eq. (16 1. In 
particular, C n $(n) = — (1/n) R n +i,o{n). 



IV. ENTANGLEMENT SPECTRUM AND 
ENTROPY MEASURES 

Once the spectral density fi(p) is known, the entropies 
can be calculated using the integral forms of Eqs. Q and 
Eq. (pH), which are 

S = - f d P fx(p) [plnp +(l-p) ln(l - p)] (17) 
Jo 

for the von Neumann entropy and 



S a — 



1 - a Jo 



dpn(p) \n[p a + (l-p) a ] 



(18) 



for the Renyi entropies. 

The spectral density can be obtained from the jump of 
hix(n,x) across the line k = ±1/2 (see, e.g., Ref. [53]): 



m(p) 



1 



d 



At: 2 p(l — p) dn 



lnx(«;, x) 



k=+(j— s)— i£ 



«=-(i-e)-i« 



. (19) 



where e is an infinitesimally small positive parameter and 
we introduced the parameterization 



A. Expansion of the spectral density fi(p) 

The spectral density may now be found by inserting 
Eq. ( 16 ) into Eq. ( 19 ). Using the symmetry of the gener- 



ating function x{~ K i x ) — X*( k t x ): we nray rewrite the 
result as 



1 



Re 







2tt 2 p(1 — p) dn 

OO |_2" J 

n=l ro= _L»j 



2inx — 2k 2 \nx + C(k) 



2imx 



K=\— i£ 



(21) 



Note that the x dependence of each term in Eq. (21) 



is known. The coefficients at nonoscillating terms are 
determined by C n ,o, so that the smooth (nonoscillating) 
part of fi(p) can be calculated as 



mO) 



1 



( ln(2x) - rip) 

lT 2 p{l-p)\ ^ ' U ' 

ERec; )0 fJ-iA x~ n 



(22) 



where r(p) is given by Eqs. (10)-(11) and the prime de- 
notes the derivative of C„,o(k) with respect to K. The 
first two terms in this expansion give the announced re- 
sult @. 
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Oscillating terms may, in turn, be collected by the 
"diagonals" C 2n ,- n +i( K ) and C 2n +i -n+lfa) witn 1 

0. 1 . '2 contributing terms of the orders x~ 2 ' and 

^,-2;-!^ reS p CC tivcly. The first two diagonals (with I = 0) 



are easy to sum. By calculating the logarithm of the se- 
ries flj-fll]) and using = 2k 3 (see Ref. [18]), we 
find 

c 2n ._„( K ) = tii^i e „[c(„-i)-c(«)] (23) 

n 

C 2n+1 ,_„( K ) = 2 l (-l) n+1 (3« 2 - 3« + 1) e »[c(«-i)-c(«)] . 

(24) 

Adding those contributions converts the continuous spec- 
trum (22 1 into a sum of delta functions. For example, 



taking into account the diagonals ( 23 1 and ( 24 1 results 



m(p) = E' 5 



$(p, X) — 7T 



d<5> 

dp ' 



(25) 



where 



$(p, x) = x+2£ ln(2x)-2( /3 (0+ f 3£ 2 - J J x^ 1 -)-© (x~ 2 ) 

(26) 

and is given by Eq. (11). 

Note that the spectrum (25) has a quasiclassical na- 



ture: the positions of quantum levels are determined by 
a quantization rule of Bohr-Sommerfeld type. The result- 
ing spectrum is regularly spaced with the average density 
given by (1/V) d&/dp. This can be explained by the fact 
that the diagonals (p3| and (24) stem, in fact, only from 



the two leading Fishcr-Hartwig branches in Eq. (13) at 
k = 1/2 (those with j = and j = —1). The spec- 
trum is thus determined from the condition that these 
two branches cancel each other, which naturally leads 
to an expression of the form ( 25 ) . Including higher- 



order Fisher-Hartwig branches produces modulations in 
the level spacing, but this effect appears only at higher 
orders in 1/x. 

Note also that this expansion breaks down close to the 
end points of the spectrum p = 1 and p = 0. In those 
regions r(p) is large and therefore the density of states 
given by Eq. (22) becomes formally negative: in fact, 
the expansion (22) is not applicable in those regions of 
p. Indeed, the expansion parameter in Eq. (22) is £/x: 
this can be seen from the (unproven) fact observed in 
Ref. [TH] that the polynomial C n ,o has degree n + 2 in k 



(and therefore in £). Thus the expansion (22) is only 
applicable at |£| <C x. Remarkably, this condition also 
guarantees the positivity of ft(p). 

Our results (|9|-(fIT]) and the quasiclassical structure of 
the spectrum are consistent with the numerical studies 
of Rcfs. [35] In particular, the smooth part of the density 
of states in the middle of the spectrum is 



£(p = 1/2) = -~ (lnx + b) + Oix- 1 ) , 
6 = ln2 -<//(()) « 2.657, 
in agreement with the findings of those works. 



(27) 



B. The entropy measures 

The von Neumann and Renyi entropies can now be 
found by substituting the expansion (21) into the inte- 



grals (17) and (18), respectively. Note that, even though 



the expansion (21) applies only at |£| <C x, we may in- 
tegrate in Eqs. (17) and (18) from £ = — oo to £ = +oo 



(corresponding to < p < 1): the contributions from 
large £ are exponentially smaller than all the terms of 
the resulting series and may be neglected. 

For the von Neumann entanglement entropy, oscillat- 
ing contributions vanish at all orders in 1 jx [the integral 
( 17 1 may be closed in the upper or lower half plane of the 
variable £, see Appendix [A]. Only nonoscillating contri- 
butions survive and may be found by replacing n(p) in 
the integral (17) by its nonoscillating part (p2|. As a 



result, we find the power series 



1 



S(x) = - ln(2x) + T + 



OO 

E 

71=1 



s n x 



(28) 



where the coefficients are given by 



d£ f TmC n J\-i^ ) m) 
cosh (7r£) 



The functions C n _o(n) may be found from the results 
reported in Ref. 18 or calculated order by order used the 
methods developed in that work. In particular, it follows 
from the results of Ref. \TE\ that C2 n +i,o( K ) are polynomi- 
als odd in k with purely imaginary coefficients. There- 
fore, all the odd terms in the expansion (28) vanish, and 
we arrive at the result (J7|. Furthermore, since C n fl{n) 
are polynomials with rational coefficients, all the coeffi- 
cients S2n are rational numbers. The first three nonzero 
coefficients can be obtained from 



~ , v 25 4 63 6 

C 40 (k) = K H K b . 

4,01 ; 16 ^ 4 



(30) 



c 6 ,o(«0 = 



889 

/ 

12 



3129 

/ 

16 



which gives the result ([8]). Following this procedure [with 
C , 2n,o('t) calculated using the method of Ref. [18], the co- 
efficients S2n may be computed to any order, one by one, 
in a straightforward way. 

The constant T is found to be 



T 



2 '■ + v 

7T 



d^ip'(0 (In [2cosh(7rO] - < taring]) 
w 0.4950179, (31) 



where the function </?(£) is defined by Eq. (11). This 
expression for T can be shown to agree with that found 
in Ref. [H 

In contrast, for the Renyi entropies, the oscillating 
parts do not vanish and can be classified in terms of the 
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poles of the integrand of Eq. (18), see Appendix |A} The 
first orders [calculated using Eqs. (23) and (24)] are 



S a {x) = 



ln(2x) + T a 



(a + l)(3a 2 -7) _ 2 



oo 



^— ' a — 1 

n,j = l 



cos(2nx) + x 



96a 3 




A V2 ^ 2a / 


1 V2 2q J. 


(2j-l) 2 " 


sin 


9 

ar 



+ o(x- 2 ). 



(32) 



These corrections reproduce the results of Ref. (in 
the corresponding continuous limit of the spin- 1/2 XX 
chain) and Refs. [TTJ A straightforward calculation of the 
constant term in Eq. ( 32 ) gives 



of the expansion to this line is however a very plausible 
conjecture, since the expansion itself is regular at this 
line and a numerical study on the more general lattice 
modeP^ indicates that, at least at the lowest orders, it 
remains valid at K = 1/2. We thus conjecture that our 
results arc in fact exact expressions for the model con- 
sidered. 
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Appendix A: Oscillating contributions to the 
entropies 



In [2 cosh(7r£a)] — a ln[2 cosh(7r£)] 
x - . (33) 

1 — a 

This expression reproduces the result found by Jin and 
Korepin in Ref. [12l 

Calculating higher-order oscillating terms in the Renyi 
entropies would require knowing higher order diagonals 
C 2 „, _„+/(«) and C 2n +i,-n+l(K) with I > 1. Although 
each of those coefficient can be separately calculated (us- 
ing the methods of Ref. [T5)) . deriving general formulas 
(valid for all n) is a tedious task, and we do not attempt 
it here. 



V. SUMMARY AND DISCUSSION 

In this paper, we have used the asymptotic expansion 
of the FCS generating function for a line segment of one- 
dimensional free fermions to determine the asymptotic 
expansion of the entanglement entropy and the entangle- 
ment spectrum in the same system. The main result is 
the asymptotic power series in 1/x for the von Neumann 
entropy. Our method also allows to construct finite-size 
corrections for the Renyi entropies (we only do it to the 
lowest order, where we reproduce the known results) and 
gives an expansion for the spectrum of the single-particle 
correlation matrix. 

Our results are based on the expansion conjectured 
(not rigorously proven) in Ref. 1181 and therefore also have 
the status of conjecture. Two elements of the proof were 
missing in Ref. 18, First, the periodicity relations on the 
expansion coefficients [which allows to convert the expan- 
sion into an explicitly periodic form ( 13 )] were not proven 
but only checked analytically up to the 15th order in 1/x. 
Second, the expansion ( 13 1— ( 14 ) was not extended to the 



line Re(ft) = 1/2: the point where the switching of the 
Fisher-Hartwig branches takes place and where we need 
the expansion for calculating the entropies. An extension 



In this appendix, we treat the oscillating (in x) terms 
in the expansions of the von Neumann and Renyi entan- 
glement entropies. For the von Neumann entropy (17), 



all the oscillating terms vanish, provided the expansion 
conjectured in Ref. [TBI is correct. For the Renyi entropies 
(18), there are oscillating terms decaying as a-dependent 



powers of x. 

Oscillating terms in the entropies are obtained by sub- 
stituting the terms of the expansion (21) with a given 



oscillation frequency m into the integrals (17) and (18). 
The integrals are further calculated by using £ as the in- 
tegration variable, integrating by parts and closing the 
integration contour in the upper (lower) half plane for 
m > (m < 0, respectively). 

In the case of the von Neumann entropy, this produces 
terms of the form 



Im 



^2,imx n— 2m 



x CVi 



cosh (7r£) 
1 
2 



4im£ In x 



(Al) 



Now the crucial ingredient of our discussion is the struc- 
ture of the coefficients C' n ^ m (n). It can be seen from 
the explicit calculation in Ref. HH] (using the Riemann- 
Hilbert method) that these coefficients have the following 
form (assuming m > 0): 

Cn,m( K ) = C n ,m('*) e ^ ^ ^ , (A2) 
Cn -m( K ) — Cn,— mi^S ^ ' ^ , (A3) 



where C(k) is defined in Eq. (15) and c„ im (/c) are some 
polynomials in k. 

From this property, it follows that, at m > 0, the co- 
efficient C n>m (l/2 — i£) has zeroes of degree 2m at all 
points £ = i(l/2 + r) for r = 0, 1, . . ., which compen- 
sate the pol es of degree two of the factor cosh~ 2 (7r£) in 
the integral Therefore the integrand is analytic in 

the upper half plane where the contour is closed, and 
the integral vanishes. Similarly, at m < 0, the coeffi- 
cient C niTO (l/2 — i£) has zeroes of degree 2m at all points 



G 



£ = —i{\/2 + r) for r = 0,1,..., the integrand is analytic 
in the lower half plane, and the integral vanishes again. 
We therefore conclude that the asymptotic expansion of 
the von Neumann entanglement entropy has the form of 
a power series in 1 jx (apart from the leading logarithm) , 
without any oscillating terms. 

In the case of the Renyi entropies, the oscillating terms 
have the form 



They contain additional poles at £ = ±(i/a)(l/2 + n). 
These poles are not compensated by zeroes of C„. m (l/2 — 
«£) and produce oscillating contributions to the entropy 
decaying as fractional (a-dependent) powers of x. A cal- 
culation of the first few terms [based on the explicit ex- 
pressions (23) and (24)] produces the result (32). 



Im 



^2imx g, — n— 2m 



r/£ 



a [tanh(-7r£) — tanh(a7r£)] 



1 - a 



x C„, m ( — — i£ ] e 



4im£ In x 



(A4) 
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